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(54) Title: FLUORESCENCE LIFETIME-BASED IMAGING AND SPECTROSCOPY IN TISSUES AND OTHER RANDOM MEDIA 
(57) Abstract 

A system and method non-invasive biomedical optical imaging 
and spectroscopy with low-level light is described. The technique 
consists of a modulated light source (120) coupled to tissue 
(100) of a patient to introduce excitation light. Fluorescent 
light emitted in response to the excitation light is detected with 
sensor (148). The AC intensity and phase of the excitation and 
detected fluorescent light is provided to a processor (160) operatively 
coupled to sensor (148). Processor (160) employs the measured 
re-emission kinetics of excitation and fluorescent light to "map" 
the spatial variation of one or more fluorescence characteristics of 
die tissue (100). The fluorescence characteristic may be provided 
by exogenous contract agents, endogenous fluorophores, or both. 
The variations is determined by solving frequency domain diffusion 
equations at a number of designated points in the tissue as part 
of a recursive estimation algorithm. Processor (160) generates 
an imaging signal in accordance with the spatial variation of the 
fluorescence characteristic for provision to an output device (164). 
The output device (164) displays an image corresponding spatial 
variation of the fluorescence characteristic which corresponds to 
tissue (100) to aid in the detection and diagnosis of disease. 




FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States party to the PCT on the front pages of pamphlets publishing international 
applications under the PCT. 



Buitina Faso 
Bulgaria 



Central African Republic 

Switzerland 
Ca»e d'lvoire 
Cameroon 

Czechoslovakia 
Czecli Republic 
Germany 
Denmark 

Spain 
Finland 



United Kingdom 



Kytgystan 

Democralic People's Republic 
of Korea 

Republic of Korea 



Liechtenstein 
Sri Lanka 
Liberia 
Lithuania 
Luxembourg 

Monaco 

Republic of Moldova 

Madagascar 

Mali 



Malawi 

Netheriands 
Norway 
New 2:ealand 

Portugal 
Romania 

Russian Federation 



Slovenia 
Slovakia 

Swaziland 
Chad 
Togo 
Tajikistan 

Trinidad and Tobago 

United States of America 

Uzbekistan 

Viet Nam 



wo 97/08538 



PCTAJS96/13658 



FLUORESCENCE LIFETIME-BASED IMAGING AND SPECTROSCOPY 
IN TISSUES AND OTHER RANDOM MEDIA 

BACKGROUND OF THE INVENTION 

The present invention relates to spectroscopic imaging of heterogeneous light 
scattering media, and more particularly, but not exclusively, relates to in vivo imaging 
of biologic tissue by mapping a fluorescence characteristic of the tissue through the 
detection of light emitted in response to excitation light from a time-varying light 
source. 

The early detection of disease promises a greater efficacy for therapeutic 
mtervention. In recent years, non-invasive techniques have been developed v/hich have 
improved the ability to provide a reliable and early diagnosis of various afflictions by 
detecting biochemical changes in the tissue of a patient. For example. Magnetic 
Resonance Imaging (MRI) has successfully monitored the relaxation of spm states of 
paramagnetic nuclei in order to provide biomedical imaging and biochemical 
spectroscopy of tissues. Unfortunately, the complexity and expense of MRI diagnostics 
limit its availability - especially as a means of routine monitoring for disease. 

Another powerful analytical technique with an increasing number of 
applications in the biological sciences is fluorescence spectroscopy. These applications 
include biomedical diagnostics, genetic sequencing, and flow cytometry. To date, there 
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are several industrial and academic institutions developing fluorescent and 
phosphorescent compounds for observing pertment metabolites and environmental 
conditions, such as Ca"*^, pH, glucose, p02, and pC02. With the development of dyes 
and photodynamic fluorescent agents which excite and re-emit in the near-infrared red 
(NIR) wavelength regime, non-invasive detection of diseased tissues located deep 
within tissues may also be possible since red excitation and re-emission light can travel 
significant distances to and from the tissue-air interface (See Wilson et al. . Time- 
Pepgndent Optical Spectroscopv and Imagine for Biomedical Ap plications. Sn 
Proceedings IEEE pp. 918-30 (1 992)). 

As exemplified by U.S. Patent Nos. 5,421,337 to Richards-Kortum et al. and 
5,452,723 to Wu et al., several investigators have suggested various procedures to 
differentiate diseased and normal tissues based on fluorescence emissions through non- 
invasive external measurements or minunally invasive endoscopic measuring 
techniques. Unfortunately, these procedures generally fail to provide a viable spatial 
imagmg procedure. One reason imaging based on fluorescence has remained elusive is 
that meaningfiil relational measurements of fluorescence characteristics from a random, 
multiply scattering media, such as tissue, are difficult to obtain. For example, 
fluorescent intensity, which is a fimction of the fluorescent compound (or fluorophore) 
concentration or "uptake," is one possible candidate for imaging; however, when this 
property is used in an optically dense medium, such as a particulate (cell) suspension, 
powder, or tissue, the local scattering and absorption properties confound measured 
fluorescent intensities. 

Besides intensity, other properties of selected fluorophores such as fluorescent 
quantum efficiency and lifetime are also sensitive to the local biochemical environment. 
As used herein, "fluorescent quantum efficiency" means the fractional number of 
fluorescent photons re-emitted for each excitation photon absorbed or the fraction of 
decay events which resuh in emission of a fluorescent photon. "Fluorescent lifetime," 
as used herein, is defined as the mean survival time of the activated fluorophore or the 
mean time between the absorption of an excitation photon and re-emission of a 
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fluorescent photon. Like intensity, measurement of these fluorescence characteristics 
is often limited to well-defined in vitro applications in the research laboratory or in 
flow cytometry where issues such as scattering, absorption, and changing fluorophore 
concentrations can be controlled or measured. Moreover, these limitations generally 
5 preclude meaningful fluorescence-based imaging of hidden tissue heterogeneities, such 
as tumors or other diseased tissue regions which cannot be detected by visual 
inspection. 

Thus, a need remains for a technique to non-invasively image multiply 
scattering tissue based on one or more fluorescence characteristics which does not 
10 require extensive information about mUinsic optical properties of the tissue, and takes 
advantage of the contrast capability offered by fluorescence yield and lifetime 
characteristics to aide in the identification of tissue heterogeneities. The present 
invention satisfies this need. 
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SUMMARY OF THE INVENTION 

The invention relates to spectroscopic imaging of heterogeneous, light 
5 scattering materials. Several aspects of the invention are novel, non-obvious, and 
provide various advantages. While the actual nature of the invention covered herein 
can only be determmed with reference to the claims appended hereto, certain features 
which are characteristic of the present invention are described briefly as follows. 

One feature of the present mvention is a technique for imaging a heterogeneous 
10 light scattering material. This process includes exposing the surface of a material to 
light from a light source and detecting an emission in response. A spatial variation of a 
fluorescence characteristic of the material is determined as a function of the emission 
with a processor. The spatial variation may be characterized by a set of values 
representative of the fluorescence characteristic as a function of position. An image is 
15 generated m accordance with the spatial variation which corresponds to the 

heterogeneous composition of the material. This technique may be applied in vivo to 
biologic tissue using external or endoscopic instrumentation to detect heterogeneities 
indicative of disease. Hie technique may include the introduction of a fluorescent 
contrast agent into the material. The fluorescence characteristic detected may be 
20 fluorescence lifetime, fluorescence quantum efficiency, a fluorophore absorption 
coefficient, fluorescent yield (a function of fluorescent quantum efficiency and 
fluorophore absorption), or another fluorescence characteristic known to those skilled 
in the art. 

In another feature of the present invention, the spatial variation of a light 
25 scattering material with a heterogeneous composition is determined by establishing an 
estimate of the optical property or fluorescence characteristic variation, determining a 
calculated emission from the material as a function of the estimate, and comparing the 
calculated emission to a detected emission to determine a corresponding error. The 
estimation of the variation is modified, the calculated emission re-determined with this 
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modified estimate, and the comparison repeated until the error reaches a desired 
minimum. An image of the material is generated from the modified estunate which 
corresponds to the heterogeneous composition. 

Accordingly, it is one object of the present invention to map a fluorescent 
5 property of a light scattering material which varies with the heterogeneous composition 
of the material to generate a corresponding image. 

It is another object of the present invention to provide a spectroscopic technique 
for non-invasively monitoring fluorescent properties of hidden tissue volumes in a 
living organism and to monitor selected metabolites of an organism in vivo. 
10 Yet another object is to provide a fluorescence imaging system and process to 

identify diseased tissue using endogenous or exogenous fluorophores as a contrast 
agent This contrast may be provided by a concentration, lifetime, or quantum 
efficiency difference between normal and diseased tissues. 

It is still another object of the present invention to provide an imaging technique 
15 and algorithm which is based on contrast of an optical property that is independent of 
local fluorophore concentration. 

Further objects, features, aspects, and advantages of the present invention will 
become apparent fi-om the drawings and description contained herein. 



- 5 - 



wo 97/08538 



PCT/US96/13658 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic illustration of a system of one embodiment of the present 
invention. 

5 FIG. 2 is a flow chart of a process utilizing the system of FIG. 1 . 

FIG. 3 is a schematic representation of a tissue phantom arrangement used to 
demonstrate various aspects of the present invention. 

FIGS. 4-7 graphically depict selected properties of equations used in the present 
invention. 

10 FIGS. 8 and 9 graphically depict convergence of sunulated determinations of 

the spatial variation of fluorescent yield and lifetime, respectively, utilizing one 
embodiment of the present invention. 

FIGS. 10-14 are images obtained from experimental examples 1-3 of the present 
invention. 

15 FIG. 15 is a schematic illustration of a system of an alternative embodunent of 

the present invention. 



- 6 - 



wo 97/08538 



PCT/US96/13658 



DESCRIPTION OF PREFERRED EMBODIMENTS 

For the purposes of promoting an understanding of the principles of the invention, 
reference will now be made to the embodiment illustrated in the drawings and specific 
5 language will be used to describe the same. It will nevertheless be understood that no 
lunitation of the scope of the mvention is thereby intended. Any alterations and further 
modifications in the described device, and any further applications of the principles of the 
invention as described herein are contemplated as would normally occur to one skilled m the 
art to which the invention relates. 

10 FIG. 1 depicts system 1 1 0 of the present invention for fluorescent imaging of 

tissue 100. Tissue 100 has surface 101 and a heterogeneous composition as represented 
by regions 102, 103 underlying surface 101. Heterogeneities 102, 103 are generally 
not detectable by visual inspection of surface 101. 

System 110 includes modulated light source 120 to supply an intensity 

15 modulated excitation light of predetermined firequency and wavelength to tissue 1 00 
via optic fiber 123. Preferably, source 120 is a laser diode of conventional design with 
a modulated output in the 1-500 MHz fi-equency range and a monochromatic ou^ut in 
the 100 to 1000 nm wavelength range. The specific wavelength is selected to excite a 
designated fluorophore in tissue 100. Beam splitter 126 may be employed to direct a 

20 small portion of the excitation signal to reference sensor 1 28 for processing purposes. 

System 1 10 also includes detection subsystem 140 which has optic fibers 143 to 
detect photons emitted from tissue 100 fi-om a number of corresponding detection sites. 
Subsystem 140 includes one or more emission sensors 148. Detection subsystem 140 
also includes an interference filter to obtain a selected emission wavelength 

25 corresponding to emission of a designated fluorophore in tissue 1 00. In one 

embodiment, subsystem 140 includes a single sensor 148 and the signals fi-om fibers 
143 are multiplexed. Preferably, sensors 128, 148 are Photo-multipUer Tubes (PMTs) 
or photodiodes but other sensor varieties, such as image intensifiers and charge-coupled 
devices, are also contemplated. 
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Sensors 128, 148 and source 120 are operatively coupled to heterodyne 
subsystem 130. Subsystem 130 is configured to obtain information about the phase, 
AC, and DC intensity of light detected with sensor 128 relative to light detected with 
the sensor 148 using conventional laser heterodyning techniques. In one embodiment, 

5 heterodyne subsystem 130 includes a signal synthesizer phase-locked to the repetition 
rate of a laser used for source 120. For this embodiment, subsystem 130 includes an 
amplifier to gain modulate sensors 128, 148 at a harmonic of a laser repetition rate 
(when a pulsed laser is used) or at the modulation fi'equency (when a modulated laser 
diode is used) plus an offset to provide the desired heterodyning. In one variation of 

10 this embodiment, an 80 MHz pulsed laser repetition rate is divided down to 1 0 MHz 
and input to the synthesizer, and a heterodyning offset of 1 00 kHz is input to the 

amplifiers for sensors 128, 148. ■, 
Sensors 128, 148 are operatively coupled to processor 160. Processor 160 

includes input/control device 162, output device 164, and memory 166. Processor 160 
15 may be an electronic circuit comprised of one or more components. Similarly, 

processor 1 60 may be comprised of digital circuitry, analog circuitry, or both. Also, 

processor 1 60 may be programmable, an integrated state machine, or a hybrid 

combination thereof Preferably, input device 162 is a keyboard or mput control of a 

conventional variety, and output device 1 66 is a Cathode Ray Tube (CRT) based video 
20 display, printer, or other image display system known to those skilled in the art. 

Memory 1 66 is preferably of the electronic (e.g. solid state), magnetic, or optical 

variety of the type readily available for use with electronic controllers or processors. 

Furthermore, Memory 166 may include an optical disk memory (CD), electromagnetic 

hard or floppy disk media, or a combination of these. j 
25 FIG. 2 depicts one mode of operation of system 1 1 0 as process 210. Process 

210 mcludes mapping the spatial variation of fluorescence yield and lifetime with 

processor 1 60 and generating an image signal m accordance with the map. Output j 
device 164 is configured to display an image in response to the image signal. Process 
210 begins by introducmg a fluorescent contrast agent into tissue 100 in stage 212. 
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This agent provides a source of fluorescent emission for detection by subsystem 240. 
The configuration of the modulated light source 120, heterodyne subsystem 130, and 
detection subsystem 140 is designed to accommodate the excitation and emission 
properties of the selected fluorescent agent In other embodiments, endogenous 
fluorophores may be alternatively or additionally employed and system 110 adjusted 
accordingly. 

In stage 214, tissue 100 is excited by light source 120 configured according to 
tiie selected fluorophore. In stage 216, tiie phase, Gobs, and log of AC intensity. Mobs, 
of the emission at each detection site "i" relative to the excitation light firom source 120 
are determined at tiie heterodyne (or offset) fi-equency. For "Di" number of detection 
sites, the detected or observed phase and AC mtensity are mdexed by "i" using tiie 
following notation: (eobs)i and (Mobs){, respectively. Processor 160 stores tiie 
relative phase and AC intensity information in memory 166. 

In stage 2 1 8, a two dimensional grid is established for an area of tissue 1 00 
selected for imaging, and a matrix of grid points is established and indexed by "j", A 
uniform seed value for the fluorescent yield, yj = ('npiax->m ^j' fluorescent 
lifetime, (T)j, at each grid point j is assigned. These values are an initial homogeneous 
guess of tiie yield and lifetime values which are modified in later stages. The term "ti " 
is the quantum efficiency of the fluorophore which varies with the environment of the 
surrounding of the fluorophore. The term " Hax->m" absorption coefficient for 
the fluorophore and is tiie product of the extinction coefficient of the fluorophore based 
on the natural log and the concentration of the fluorophore. As a result, the yield, y = 
Tll^ax->.m' influenced by the surrounding metabolism and the uptake of the 
fluorophore. The uptake of certain known fluorophores vary witii the type and 
condition of host tissue, providing anotiier fluorescence characteristic useful to detect 
disease. The contrast provided by these properties is largely independent of 
fluorophore concentration. The initial estimate of fluorescent yield and lifetime are 
stored in memory 166 by processor 160 for later use. 
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After establishing this initial estimate of the fluorescence characteristics of 
yield, Wax-^m' lifetime, t, processing loop 220 is entered in stage 230. 
Preferably, the stages of processing loop 220 are executed by processor 160 via 
preprogrammed software, dedicated hardware, or a combination of both as appropriate. 
To aid in understanding various mathematical aspects of process 210 and loop 220, the 
following table of selected variables is listed: 



c velocity of light; 

D(r) optical diffusion-coefficient; 

Di number of detection sites; 

f modulation frequency; 

I identity matrix; 

i detection site index; 

J Jacobian matrix relating the sensitivity at each grid point, j, to 

the response at each detection site; 

j grid point index; 

Jj i individual elements of the Jacobian matrix J; 

k source index; 

M log of AC intensity of modulated fluorescent light position; 

m index to multiple modulation frequencies; 

n average index of refraction; 

r position (in two or three dimensions); 

Sk number of modulated light sources; 

S(r, ffl) source term for the modulated light at position r and frequency 
o; ■ 

Greek 

merit function representing the least squares error; 
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complex number representing photon flux in the frequency 
domain at position r and frequency co; 
quantum efficiency of fluorescent probe or dye; 
average absorption coefficient; 

absorption coefficient of the fluorescence light by both the non- 
fluorescing chromophores and fluorophore; 
absorption coefficient of the excitation light by both the non- 
fluorescing chromophores and fluorophore; 
adsoiption coefficient due to non-fluorescing chromophores; 
adsorption coefficient of excitation light by fluorophores; 
effective scatting coefficient; 
phase-shift of one modulated light wave to another; 
lifetime of activated probe or dye at location r; 
angular modulation frequency, given by 2nf; 

observed or experimental data; 
excitation light; and 
fluorescence or emission light. 
20 

In stage 230, phase and relative AC intensity at each detection site "i" is 
calculated as a function of the initial estimates of yield and lifetime for each grid point 
j. The calculated phase and mtensity are represented at each detection site i as (9jn)i 
and (Min)i, respectively. The values for (0in)i and (Mm)! are determined using the 
25 diffusion equation approximation of the radiative transport equation. The diffusion 
equation approximation describes the spatial and temporal transport of light in tissues 
or multiply scattering media. A coupled frequency domam dififiision equation can be 
used to predict the excitation and emission fluence rates, O^Cr, ©) and (r, oo), 
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respectively, at any location r within the selected grid of tissue 100 via equations (1) 
and (2): 

V . [Dx(r)V<Dx(r, o>)]-[nax(r)+i CD/Cn]a>x(r, ©) +Sx(r,(D) = 0 (1) 

5 

V . [Dm(r)V(I)m(r, a)]-lHaJiT)+i (o/Cj,]Om(r, co) +Sin(r,o>) = 0 (2) 

The source teim for the excitation light Sx(r,ffl) is due to the sinusoidally modulated 
10 light at an angular frequency © = 2(nf where f is typically in the MHz frequency range. 
The first term in both of the diffusion equations (1) and (2) represents the diffusive or 
"random-walk" transport of light where Dx^m is the optical diffiision coefficient of 
equation (3) as follows: 

15 Dx,m=[3(Hax,m+^^'sx,m)]■^ (3) 

and and n'g are the absorption and isotropic scattering coefficients, respectively, for 
tissue 100, the medium of interest. The optical properties are dependent on the 
wavelength of light and thus are different for the excitation light from source 120 

20 (subscript x) and fluorescent emission detected with subsystem 140 (subscript m). The 
total absorption coefficient at the excitation wavelength, ^ax» ^o contributions 
from non-fluorescing chromophores as well as from fluorophores responsive to the 
excitation wavelength. The total absorption coefficient is given by the sum of 
absorption coefficients due to non-fluorescmg chromophores, Hax->c' fluorophores 

25 ^ax->m* Generally it may be assumed that the absorption experienced at the 

fluorescent wavelength is due primarily to non-fluorescing chromophores. The velocity 
of light in tissue is Cj^=dn where n is the average index of refraction. The source term 
for the fluorescent emission is dependent on the excitation light fluence, <J>^(t, ca) and 
is given by equation (4) as follows: 

- 12 - 
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SmM = ilMax-^m W*^x(r, - icoT(r))/(l + a)2T(r)2)] 



(4) 



This terai arises from the Fourier transform of the fluorescence decay term in the time 
5 domain following an incident pulse of excitation light where: x is the fluorophore 



product of the extinction coefficient based on natural log and the concentration of the 
fluorophore in the ground state. As previously indicated, the combined product, 
^max-^m' t^rm^ fluorescent yield, y, and is proportional to the generated 

10 fluorescence fluence. Substitution of equation (4) into equation (2) facilitates 

determination of (P^i fo"" ^^^^ Srid point "j." The solution of the diffiision equations 
(1) and (2) for the two-dimensional area defined by the grid points "j" may be readily 
extended to three dunensions to estimate spatial variation of one or more fluorescence 
characteristics in a selected volume with "r" corresponding to position in three 

15 dimensions. 

Both diffiision equations (1) and (2) are linear complex elliptic equations that 
can be solved as boundary value problems for the complex quantities <Dx(r, o) and 
Oni(r, ©). This solution employs the method of finite differences to create 
corresponding finite difference equations. These difference equations are utilized to 
20 obtain an approximate solution at each grid point, j. This method of solution is 

described in other contexts in Fulton et al., Multigrid Met hod for Elliptic Problems. A 
Review. 1 14 American Meteorological Society pp. 943-59 (May 1986); and B.W. 
Pogue et al., Initial Assessment of a Simple System for Frequencv Domain Diffiise 



25 preferred method of performing this solution is with the MUDPACK routines described 
in Adams, J.C., MTTDPACK: Multi^d Por table Fortran Software for the Efficient 
Solution of Lin ear Elliptic Partial Differential Fxiuations. 34 App. Math Comp. p. 133 
(1989). For the solution of the diffiision equations, it is assumed that <I'ni,x(r.<o) = 0 on 
the surface 101 of tissue 100 which is known as the zero fluence boundary condition. It 



lifetime, r\ is the quantum efficiency, and the absorption coefficient, ^ax->m' 




40 Physics in Medicine and Biology pp.1709-1729 (1995). One 
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should be recognized that other bouridaiy conditions may be selected and the method of 
solution varied accordingly. 

The difiusion equations (1) and (2) may be solved for a complex number for 
at each grid point, j. The detected signal at the surface is proportional to the 
normal component of the gradient of the photon fluence. To approximate the signal at 
detector site "i" located on surface 101 of tissue 100, the value at an intemal grid 
point closest to the site is selected which follows from the relationship that the normal 
component of the photon fluence gradient is proportional to just inside surface 101. 
The calculated phase-lag, ^nd the log of AC intensity, Mm* at the detection sites 
"Di" are calculated from the imagmary and real parts of the complex value with 
respect to the phase and the AC intensity of source 120. 

The diffusion equations (1) and (2) provide insight into the sensitivity of 
changing the fluorescent optical properties of tissue 100 on Gm and Mm measured at 
the detector sites i. This insight results from a series of calculations while fixing 
various parameters of the diffusion equations (1) and (2). These calculations assume 
circular tissue phantom 300 with an embedded, heterogeneity 302 hidden in phantom 
background 303 as illustrated in FIG. 3, A two-dimensional grid is established for 
phantom 300 and may easily be expanded to three dimensions. Under these simulated 
conditions, a large value is assigned to absorption coefficients for both excitation and 
fluorescent light at all grid points outside the simulated tissue phantom. The four 
sources S1-S4 of FIG. 3 (Sk = 4) are simulated by assigning an arbitrary complex 
number at a grid pomt near the surface closest to each source. The twenty detection 
sites D1-D20 of FIG. 3 (Di = 20) are simulated by using the calculated values 
determined from at the grid point "j" closest to the detection site. The simulated 
solutions to diffusion equations (1) and (2) were obtained in two dimensions for a 65 x 
65 grid covering a 100mm diameter circular tissue phantom 300 with a circular, 
embedded heterogeneity having a 30mm diameter and located at the center of the tissue 
phantom 300 (this location differs slightly from the configuration of heterogeneity 302 
of FIG. 3). The simulated measurements of fluorescent phase-shift and AC mtensity 
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are reported for 20, equally spaced, circumferentially located detection sites D1-D20. 
The modulation frequency, f, was set equal to 150 MHz. The optical properties of the 
heterogeneity and the background are shown in Table 1 as follows: 
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frequency 


(mm"l) 


(nmr^) 
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background 


background 
(ns) 


(MHz) 


^^ax-»c + 


0.0 


1.0 


0.0 


1.0x10-5 


1.0 


150.0 

















In order to evaluate the influence of imax-^m' ®m computed at 

each detection site D1-D20 as the value of ^max-^m ^ *® heterogeneity increased 
from 10"^ mm"^ to 10"^ mm"^ and as TiHax-^m ^ background 303 was maintained 

15 constant. The lifetime, x, was set equal to 1 ns for both the object and the background 
causing conti^ due to differences in "ntiax-^m* ^® °f ^m ^m ^re shown in 
FIGS. 4 and 5 respectively for one active source SI. As T^ax-^m of heterogeneity 102 
increases to higher values, the AC intensity approaches an upper limit similar to what is 
expected in dilute non-scattering solutions. FIG. 5 shows how the fluorescent phase- 

20 shift, Gm, decreases as the absorption coefficient due to the fluorophore, Hax->m 

decreased 10 to 100 times die background. From these sunulations, M^ appears to be 
directiy dependent upon changes in TiHax->m °^ ^ simulated tissue heterogeneity 102 
whereas Om is indirectiy dependent on Tmax->m ^° changes in photon migration. 
In order to evaluate the influence of t, Gm and M^ were calculated at each 

25 detection site D1-D20 as tiie values of x in the heterogeneity varied from 10"^ ns to 1 0^ 
ns and the value of t in the background was held at 1 ns. The background imax->m 
was set to lO'^mm"^ and TiHax->m *® heterogeneity was set to lO-^mm"^. As 
shown in Fig. 6, the detected AC intensity increases as x decreases. Fig. 7 illustrates 
die values of the fluorescent phase-shift at each detection site as the lifetime of tiie 
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heterogeneity is changed from 0. 1 ns to 1000 ns. At a given modulation frequency (150 
MHz in this calculation), 6^ first inCTeases, reaches a maximum and then subsequently 
decreases as t is increased from 0.1 ns to 1000 ns. Therefore, both Gm and at each 
detection site D1-D20 appear to be directly influenced by the value of lifetime in the 
heterogeneity. 

Referring back to FIG. 2, in stage 240, the calculated emission phase and 
intensity, (9ni)i and (Mm)!, are compared to the measured emission phase and intensity, 
(0obs)i and (Mobs)i. for each detection site "i" to identify a difference or "error" 
between the measured and calculated values. Because (^IHax-^m^j '"^P^'^^ (j^m)i » 
comparison is posed in the form of the merit function Xji^ of equation (5) as follows: 

Sk Di 

Xjx2 = (l/Sk) 2 (1/Di) S [((Mobs)i-(Mi„)iyoM]2 (5) 
k=l i=l 

where oj^ is the typical standard deviation of noise in M^, taken to be 0,01; Sk = 
number of excitation source sites indexed by k; and Di = number of detection sites 
indexed by i. The goal of the algorithm is to minimize %^ by appropriate updates of 
(^l^ax-» w)'y ^^^^ ^ iidtial update of (imax-> m^j ^°*^er merit function in terms of 
(T)j participates in the comparison of stage 240. This merit function, Xx^ , is presented 
as equation (6) as follows: 

Sk Di 

Xx2 = (l/Sk) Z (1/Di) S [((Mobs)i-(Mni)i)/aM]2 + [((eobs)i-(0m)iy^e]2 (6) 
k=l i=l 

where oq is the typical standard deviation of noise in (0ni)i, taken to be 1 degree; Sk = 
number of excitation source sites indexed to k; and Di = number of detection sites 
indexed to i. Since the lifetime influences both (6jn)i and (Mjn)i, the phase and AC 
intensity are used in equation (6). 
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After the comparison of stage 240 is performed by calculating the merit 
functions i^, %^ , control flows to conditional 250 to test whether the comparison of 
the observed values, (Gobs)! (Mobs)i> to calculated values (9jji)i and (Mj^)} 
via the merit functions meets a selected convergence criteria. This criteria corresponds 
to the degree of tolerable error in determining the yield and lifetime values. For one 
embodiment, convergence is achieved when any of the following three quantities, (i) 
X^, (ii) change in in successive iterations of loop 220, (ii) relative change in x^ in 
successive iterations of loop 220 is less than a predetermined threshold value of 1.0 x 
10"2, In other embodunents a different comparison calculation and associated 
conditional may be employed as would occur to one skilled in the art. If conditional 
250 is satisfied, control flows to stage 270 and loop 220 is exited; however, if the 
criteria is not satisfied, execution of loop 220 continues in stage 260. 

In stage 260, the yield, (y)j = (TlHax->m>'j' lifetime, (T)j, for each grid point 
j is updated so that these values may reach the mmimum error correspondmg to the 
comparison stage 240 and conditional 250 test. In order to update these values, 
Jacobian matrices are used which describe the sensitivity of the response at each 
detection position i to changes m (y)j = (il^iax->m^j' lifetime, (T)j, at each grid 
pomt,j. Three Jacobian matrices are employed: J (M, Tm^x-). m^' J(M,x);and 1(0, 
t). The elements Ti j of these Jacobian matricies are given by Ji j = 
[aMi/(a(Tmax_>ni)jl' -^^J " [SMj/drj]; and Jij = [ae/axj], respectively. These elements 
may be calculated by solving the diffusion (1) and (2) four tunes for each grid point, j 
to obtain U^-^ and calculated with (x )j and (x+ 6x)j and vwth (tmax-^m^j 
(Tmax->m "'" ^^f^x->m)j- squares minimization, the update to yield and 

lifetime is calculated. In one preferred embodiment, this updatmg algorithm is adapted 
- 17 - 
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from an algorithm used to reconstruct images obtained by electrical impedance 
tomography like the algorithm suggested by Yorkey, et al. Comparing reconstruction 
Algorithms for Electrical Impedance Tomography. 34 Transactions in Biomedical 
Engineering pp. 843-52 (1 987). The Jacobian Matrices are used to solve for update 
5 vectors, [A"tmajj_^jjj ] and [Ax], to estimated yield and lifetime vectors, [Ti|iax-> m ^ 
and [t], respectively. These vectors are of a dimension corresponding to the number of 
grid points. At each iteration through loop 220, the following Jacobian equations (7) 
and (8) are solved to determine the update for the estimated yield and lifetime vectors: 



J(M.m )'J(M.m ) 



J(M.m )^ 



-(Mmota - Mm) 



(7) 



j(m.t)'''j(m.t) j(e.T)'''j(e.T) , • 



^2 obs ^2 obs 



(8) 



20 ^rn^bs 2nd Mjji are the observed and calculated vectors of the log of AC intensity at 
each of the i detection sites, respectively, Q^^^ and 9^ are the observed and 
calculated vectors of the phase lag at each of the i detection sites, respectively. Due to 
the ill-conditioned nature of the Jacobian matrices, the terms X^I or X'2l are added as 
part of a Marquardt minimization scheme where I is an identity matrix. The parameters 

25 X\ or ^2 are adjusted via a Maquardt-Levenberg type algorithm of the type disclosed in 
Press et al., Numerical Re cipes: The Art of Scientific Computing. (Cambridge 
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University Press, 1992). Conventional numerical methods are employed to solve the 
simultaneous linear algebraic equations resulting from the Jacobian matrix equations 
(7) and (8). The Jacobian matrices are re-calculated at each iteration through loop 220. 
It has been foimd that equations (7) and (8) provide a way to select appropriate changes 
to the yield and lifetime estimates; however, other numerical approaches to recursively 
iterate to acceptable estunates as would occur to one skilled in the art are also 
contemplated. Once the update is complete, control returns to stage 230. 

If the convergence criteria is satisfied in conditional 250, then estimation of 
yield and lifetime for the grid points has reached an acceptable minimum and control 
flows to stage 270. In stage 270 an image signal is generated by processor 160 from the 
spatial variation of the yield and/or lifetime fluorescence characteristics. This image 
signal is sent to output device 164 which displays an image in response. Because the 
fluorescence characteristics of yield and lifetime typically vary with the biologic 
environment of the fluorophore, this image is generally indicative of tissue variation 
and offers the capability to detect heterogeneities 102, 103. . For example, laser diodes 
capable of supplying Near infrared (NIR) light that can penetrate tissue several 
centimeters, and fluorescent contrast agents responsive to NIR light may be used to 
provide a viable imaging system. In one embodiment, this system is adapted for use 
with an endoscope. 

Besides yield and lifetime, the spatial variation of other fluorescence 
characteristics usefiil to distinguish diseased tissues may be mapped using tiie diffusion 
equations (1) and (2). Such alternative fluorescence characteristics include, but are not 
limited to, quantum efficiency r\ and/or fluorescent absorption coefficient Rax->m 
determined as separate properties independent of the yield product. 

In another embodiment of the present invention, the photon fluence equation 
and Jacobian estimation process is adapted to determine a map of a designated 
fluorophore uptake concentration. For this embodiment, a first map of chromophore 
adsorption coefficients piax-*c scattering coefficients are determined in the 
absence of the designated fluorophore by estimating the chromophore adsorption 
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coefficient ^ax-n; scattering coefScient \i\ at each grid point j in place of the 
yield and lifetime estimates. Dif&ision equation (1) for Ox(r, co) may be used in 
conjunction mih modified Jacobian equations (7) and (8) to create this first map. The 
modification substitutes the chromophore adsorption and scattering coefficients in place 
of the yield and after adsqjtation to accommodate these new characteristics as follows: 



J(M_,fi. )^J(M„,n. ) J(e_.|x- )Tj(e„ 



of 



(9) 



J2 „2 2 



(10) 



The elements of the four Jacobian matrices employed, J(Mx,Ma,_ ). J(Mx,tis). J(6x,Ha ). 



and J(ex,Hs) sre given by jys 



respectively. Updates to the absorption and scattering map were conducted to minimize the 
merit function x^: 

. ^. . „ .^2 o ^2 

= 



.A°f A|f^xobs,i-Mx,ir re,,,3,-e,,^^ 



v^ere n. = Sk and n^ = Di. 
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After generating the first map, the designated fluorescent contrast agent is 
introduced, and the total adsorption coefficient li^^ is determined by substituting ji^x 
in place of |iax->c equations (9)-(l 1) to obtain a second map of the total adsorption 
coefficient. Noting that " ^^ax-^m ^^x-^ c *® uptake of the 

5 fluorescing contrast agent is directiy proportional to fiax-j-m' ^P^^ concentration may 
be mapped by determining a difference between the adsorption coefficient variations 
for the first and second maps. This "difference map" may then be used to generate an 
image corresponding to the uptake concentration. 

Another alternative embodiment measures the emission responsive to each of a 
10 number of light source modulation firequenciesf The total number of different 

firequencies employed is designated Mf. To obtain this additional data, an iteration of 

loop 220 is performed for each fi-equency f indexed to m. The number of sources, Sk | 
and detection sites Di are indexed to k and i, respectively. This additional data may be 
used to enhance imaging results obtained with system 1 1 0 or to permit reduction of the 
15 number of detection sites or excitation source sites in the evaluation. A representative 
merit function corresponding to this additional data is given in equation (12) as follows: 

20 

Mf Sk Di 

(IMf) 2 (1/Sk) 2 (1/Di) 2 [((Mobs)i - (Mm)iymP- + [((0obs)i " (0m)i)/<^0]2 (12) 
m=l k=l 1=1 

25 Besides fluorescence yield and lifetime, the multi-frequency method can be employed 
to map other optical characteristics of interest. Besides a sinusoidally modulated light 
source, the present invention may be adapted to operate with a pulsed or other time- 
varying excitation light source in alternative embodiments. 

FIG. 15 depicts an optical system 410 of another embodiment of the present 

30 invention. This system includes modulated light source 420 wth laser driver 422, 
operatively coupled laser diode 424, and reference frequency generator 426. Source 
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420 is configured to deliver modulated light to tissue phantom 400, and the re-emitted 
light from the phantom is focused onto a gain modulated image intensifier 430 via 
50mm lens 432. Intensifier 430 includes a photocathode face which converts photons 
to electrons, a Multi-Channel Plate (MCP) which multiplies the electronic signal by 

5 avalanche multiplication, and a phosphorescent screen which converts electrons into an 
optical image. Preferably, intensifier 430 is a fast intensifier, of the variety 
manufactured by Litton Electronics, Inc., which enables modulation by applying a DC 
bias and an RF signal from amplifier 428 between the photocathode and the MCP. For 
this example, the modulation of the image from intensifier 430 is phase-locked to the 

10 laser diode 424 by a 1 0 MHz output signal from synthesizer 426. By modulating the 
laser diode 424 and the image intensifier 430 at the same frequency, a steady-state 
image results on the phosphor screen. U.S. Patent No. 5,213,105 to Gratton et al. 
provides additional background concerning certain aspects of this technique. The 
uhage from the phosphor screen is focused through interference filter 433 on a Charge 

15 Coupled Device (CCD) camera 434 via 1 50 mm macro lens 436. Camera 434 has a 
512x512 array of CCD detectors configured to provide a corresponding pixelated 
image. Camera 434 is operatively coupled to processor 460 of a similar configuration 
to processor 160 previously described. 

Following each acquired image, a phase delay between the image intensifier 430 

20 and the laser diode 424 is induced by stepping the phase of the image intensifier 430 to 
values between 0 and 360 degrees with the frequency synthesizer 452 under the control 
of processor 460. Since the gain modulation of image intensifier 430 and laser diode 
424 occurs at the same frequency, homodyning results in a steady phosphorescent 
image on intensifier 430 which is dependent upon phase. Preferably, control between 

25 synthesizer 452 and processor 460 is obtained by a conventional GPIB interface. 

Images from the phosphorescent screen of the image mtensifier 430 are then gathered at 
each phase delay. The incremental phase delayed images are then used to generate a 
map of phase-shift and intensity modulation ratio between the excitation and re-emitted 
light from phantom 400. By applying interference or appropriate optical 
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filters, the emission light may be selectively separated from the excitation light and 
measured. Camera 434 output may be processed by processor 460 using process 210. 

The present invention will be further described with reference to the following 
specific examples 1-3. It will be understood that these examples are illustrative and not 
restrictive in nature. Examples 1-3 involve the computer simulation of the process 
210. Simulations of this kind, includmg the simulation of tissue, are an acceptable 
means of demonstrating fluorescent spectroscopic imaging performance to those skilled 
in the art. The examples use simulated values obtained by solving the diffiision 
equations (I) and (2 ) for Gm and Mjjj under the conditions of table 2 as follows: 



Case 






Rs, or 


1 




Gaussian 


Gaussian 










(background) 


(background) 


Noise 
in log of AC 


Noise 
in phase 




(mm-1) 


(mm-1) 


(mm-1) 


(ns) 


(mm-i) 


intensity 
OM 


CQ (degrees) 


5.1 


0.0 


0.0 


1.0 


10.0 


1.0x10-5 


0.01 


0.1 


5.2 


1.0 X 10-3 


0.0 


1.0 


10.0 


1.0 X 10-5 


0.01 


0.1 


5.3 


0.0 


0-0 


1.0 


10.0 


1.0 X 10-5 


0.01 


1.0 



Table 2 

The examples simulate tissue phantom 300 of FIG. 3 having a 100mm diameter. 
Values of 9^^ and Mjjj were computed at each of the D1-D20 detection sites of FIG. 3 
m response to the 4 modulated light sources S1-S4 located at the periphery. The 
excitation light modulation frequency f was simulated at 150 MHz. DifiEiision 
equations (1) and (2) were solved to provide 80 simulated values of and Mj^ 
corresponding to the various combinations of detection and source sites (Sk * Di = 
4x20 = 80). Gaussian noise with a standard deviation of 0.1 degrees (or a liberal 1 
degree) in and 1% in Mjn were superimposed on the diffusion equation solutions. 
Adapted MUDPACK routines were used to solve the diffusion equations (1) and (2) on 
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a SunSparclO computer. These obtained data sets were used as simulated input data to 
process 210 for examples 1-3. The results are shown in tables 3 and 4 are as follows: 



5 


Case 


Area, object I 
(mm2) 


Location, object 1 
(x,y), (mm, mm) 


Area, object 2 
(mm2) 


Location, object 2 
(x,y), (mm, mm) 




5.1 


706.0 (expected) 
742.2 (obtained) 


(60,60) (expected) 
(60.8,58.5) (obtained) 


not applicable 


not applicable 


10 


5.2 


706.0 (expected) 

703.1 (obtained) 


(60,60) (expected) 
(59,4,58.3) (obtained) 


not applicable 


not applicable 




5.3 


314.1 (expected) 
38 1.0 (obtained) 


(32.3,67.7) (expected) 
(34.0,67.7) (obtained) 


314.1 (expected) 
342.0 (obtained) 


(67.7,32.3) (expected) 
(65.0,35.0) (obtained) 



Tables 

15 



Case 


^{^a,^„ (object) 


T (object) 




(mm-1) 


(ns) 


5.1 


1.0 X 10-3 (expected) 


1.0 (expected) 




0.93 X 10-3 (obtained) 


1.03 (obtained) 


5.2 


1.0 X 10-3 (expected) 


1.0 (expected) 




0.8 X 10-3 (obtained) 


0.7 (obtained) . 


5.3 


(top left object): 


(top left object): 




1.0 X 10-3 (expected) 


1.0 (expected) 




2 X 10-3 (obtained) 


4.1 (obtained) 




(bottom right object): 


(bottom right object): 




2.0 X 10-3 (expected) 


2,0 (expected) 




1.8 X 10-3 (obtained) 


3.5 (obtained) 



Table 4 



wo 97/08538 



PCT/US96/13658 



Example 1 

Example 1 reconstructs fluorescent yield and lifetime with no absorption due to 
non-fluorescing chromophores. To simulate the experimental data for this example, 
5 the fluorescent yield, (^l^ax-^m^j' *® background and the heterogeneity 302 were 
chosen as 1 x 10'^ mm"^ and 1x10"^ mm'^ respectively and the fluorescence lifetime, 
(x )j, for the background and the heterogeneity 302 chosen as 10 ns and 1 ns 
respectively. During the execution of loop 220, no a priori knowledge of either the 
heterogeneity 302 location or the backgroimd fluorescence properties was assumed and 

10 a uniform guess of Ix 10"^ mm"^ and 10 ns was given for the fluorescence yield, 

('HKax-j.m^j' lifetime, (x )j, respectively. Convergence was achieved in less than 50 
iterations of Loop 220 (computational time on a SunSparclO: 2 hours) for a two 
dimensional 17x17 grid. The average values of Wa.x-^m *® ^rid points 

which occupy the simulated object converge within 50 iterations to m^ax-i-m ~ ^"^^ ^ 

15 10"3 rnm-l and x =1.03 ns are illustrated in FIGS. 8 and 9, respectively. FIGS. 10 and 
1 1 illustrate the reconstructed images from the mapped values of Wax-^m 
respectively, and are representative of the expected images. The images were smoothed 
by interpolation in examples 1-3 to remove spurious points which had unphysically 
high values, but were surround by values within a physically achievable range. These 

20 spurious values were replaced by the average background fluorescence yield and 
lifetime obtained from simulation of loop 220. 

The average values of "nUax-^ m *® points which occupy the simulated 
background converge within 50 iterations to 9 x 10"^ mm"l. The value of the 
background converges to 5.4 ns. The dependence of the final hnages on the choice of 

25 the initial guess was examined by providing an initial uniform guess of 1 x 10"^ mm'^ 
and 10ns for (^^ax->m)j' lifetime, (x)j, respectively. This resulted in shnilar 
images to those obtained in FIGS. 10 and 1 1 . 

The location of heterogeneity 302 was identified as consisting of all the grid 
points with Tmax_> m ^S)^^^ (arbitrarily chosen) of the peak value of the 
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^max-^rn (FIG- 10)- The average of the coordinates of all the identified object grid 
points was the position (60.8, 58.5) which is close to position (60, 60) that was used to 
simulate the experimental data. As listed in Table 3, the area of the heterogeneity based 
upon our arbitrary definition for identification was 72 nun^, close to that used to 
generate our simulated experimental data. 

Example 2 

Example 2 reconstructs fluorescent yield and lifetime with a simulated 
chromophore absorption configured to mimic tissue. The same hidden heterogeneity as 
well as optical parameters and simulation equipment were used as described in 
Example 1 except that a imiform background chromophore absorption coefficient, 
Max-)- °^ ^ ^^^'^ *° generate the simulated experimental data. While 

excitation light propagation was not employed for image reconstruction, we considered 
this optical property known to estimate the best possible performance for inverse image 
reconstruction under physiological conditions. The two-dimensional reconstructed 
spatial map ofthe fluorescence yield, (Tmax_).in)j' and lifetime, (T)j, are shown in 
FIGS. 12 and 13, respectively. As shown in Table 3, the mean value of location of the 
object according to our criterion based on imax_>m occurred as position (59.4, 58.3) 
consistent with the conditions used to simulate the experimental data. The dimension 
of the heterogeneity based upon our arbitrary definition for identification (all grid 
points with T]liax-^m ^8^^^^ *an 35% of the maximum) were 703 mm^ which is close 
to that used to generate our simulated experimental data. The average values of 
Tl^ax-^m ""^ *® S"*^ points which occupy the simulated object converge within 
50 iterations to the values of Tifiax-^ni ~ ^-^ ^^^'^ and x = 0.7 ns consistent 
with the values used to generate the simulated experimental data (see Table 3). The 
average values of '»max-»m points which occupy the simulated 

background converge within 50 iterations to values similar to that reported for Example 
1. 
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Example 3 

Example 3 simulated two hidden heterogeneities in the tissue phantom (not 
5 shown in FIG. 3). In this case, the same optical parameters were used as described in 
example 1 except that the fluorescence yield imax-j-m objects 1 and 2 was 
chosen as 1 x lO'^mm"^ and 2 x lO'^mm"^ respectively and lifetime t for the 
heterogeneities chosen as 1 ns and 2 ns, respectively. A 33 x33 grid was employed 
instead of a 17 xl 7 grid. An image corresponding to the mapping of yield is depicted 
10 in FIG. 14. 

All publications and patent applications cited in this specification are herein 
incorporated by reference as if each individual publication or patent application were 
specifically and individually mdicated to be incorporated by reference. While the 
15 invention has been illustrated and described in detail in die drawings and foregoing 

description, the same is to be considered as illustrative and not restrictive in character, it 
being understood that only the preferred embodiment has been shown and described 
and that all changes and modifications that come within the spirit of the invention are 
desired to be protected. 
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CLAIMS 

What is claimed is: 

1 . A method of imaging, comprising the steps of: 

(a) exposing a surface of a light scattering material to an excitation light from a light 
source, the material having a heterogeneous composition under the surface; 

(b) detecting a fluorescent emission from the material in response to step (a); 

(c) establishing an estimate of spatial variation of a fluorescence characteristic of the 
material; 

(d) determining a calculated emission as a function of the estimate; 

(e) comparing the calculated emission to the emission detected in step (b) to determine 
an error; 

(f) providing a modified estimate of spatial variation of the fluorescence characteristic 
and repeating steps (d) through (f) xmtil the error reaches a desired minimum; and 

(g) generating an image of the material from the modified estimate, the image 
corresponding to the heterogeneous composition of the material. 

2. The method of claim 1, further comprising introducing a fluorescent contrast agent into 
the material. 

3 . The method of claim 1 , wherein step (f) includes applying a Jacobian matrix. 

4. The method of claim 1 , wherem the expected emission is determined as a ftmction of a 
difiiision equation. 

5. The method of claim 1, wherein step (e) includes comparing intensity and phase of the 
calculated emission to intensity and phase of the emission detected in step (b). 
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6. The method of claim 1, wherein the material includes a fluorescent contrast agent and 
the fluorescence characteristic is a function of at least one of fluorescence quantum 
efBciency, fluorescence lifetime, and concentration of the fluorescent contrast agent. 

7. The method of claim 6, wherein the fluorescence characteristic is a function of 
fluorescence quantum efBciency, the light source is intensity modulated at a predetermined 
frequency, step (e) includes comparing AC mtensity and phase of the calculated emission to 
intensity and phase of the emission detected in step (b), step (f) includes applying a Jacobian 
operator, and the expected emission is determined as a function of a photon fluence rate 
within the tissue. 

8. A method of imaging, comprising: 

(a) exposing a living biologic tissue to an excitation light from a light source; 

(b) detecting a fluorescent emission from the tissue in response to said exposing; 

(c) determining a spatial variation of a fluorescence characteristic within the tissue as a 
function of the emission with a processor; and 

(d) generating an image of the tissue in accordance with the spatial variation. 

9. The method of claim 8, further comprismg introducing a fluorophore into the tissue. 

10. The method of claim 8, wherein the fluorescence characteristic is a function of at least 
one of fluorescence lifetime, fluorescence quantum efficiency, or fluorescence absorption. 

1 1 . The method of claim 9, wherein said exposing includes positioning the source adjacent 
a surface of the tissue, and said detecting includes sensing emitted light at the surface. 

12. The method of claim 8, wherein said ejqjosing includes positioning a number of 
modulated light sources adjacent a sur^ce of the tissue and said detecting includes sensing 
emitted light at a number of locations along the surface. 
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13. The method of claim 8, wherein said detemiining includes (i) establishing an estimate 
of spatial variation of the fluorescence characteristic, (ii) determming a calculated emission as 
a ftmction of the estimate, (iii) comparing the calculated emission to the fluorescent emission 
of said detecting to determine an error, (iv) providing a modified estimate of spatial variation 
of the fluorescence characteristic, and (v) repeating said determining and said comparing with 
the modified estimate until the error reaches a desired minimum. 

14. A method of imaging, comprising: 

(a) introducing a fluorescent agent into a light scattering biologic tissue with a 
heterogeneous composition; 

(b) exposing a surface of the tissue to light fi-om a light source to excite the agent; 

(c) detecting an emission of light firom the tissue in response to said exposing; 

(d) determining a number of values representative of a fluorescence characteristic of 
the tissue as a fimction of position and the emission; and 

(e) generating an image in accordance with the values, the image corresponding to the 
heterogeneous composition of the tissue. 

15. The method of claim 14, wherein the fluorescence characteristic is a function of 
fluorescence quantum efficiency. 

1 6. The method of claim 1 4, wherein the fluorescence characteristic corresponds to 
fluorescence lifetime. 

1 7. The method of claim 14, wherein said detecting includes sensing fluorescent light 
emitted in response to said exposing at a number of locations along the surface. 

1 8 . The method of claim 1 4, wherein said exposing includes providing light fi-om a number 
of sources spaced along the surface. 
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19. The method of claim 1 4, wherein said exposing includes exciting the fluorescent agent 
in the tissue with modulated light at a number of different frequencies, and the values are 
determined as a function of the different frequencies. 

20. The method of claim 1 4, wherein said determining includes (i) establishing an estimate 
of the values, (ii) determinmg a calculated emission as a fimction of the estimate, (iii) 
comparing the calculated emission to the emission of said detecting to determine an error, (iv) 
providmg a modified estimate of tiie values, and (v) repeating said determming and said 
comparing with the modified estimate until tiie error reaches a desired minimum. 

21. A system for imaging a light scattering tissue having a heterogeneous composition and 
containing a fluorophore, comprising: 

(a) a light source adapted to excite the fluorophore; 

(b) a sensor configured to provide a detected light signal corresponding to a fluorescent 
emission from the tissue in response to light from said source; 

(c) a processor operatively coupled to said sensor and responsive to said detected light 
signal to provide a number of values representative of a fluorescence characteristic of tiie 
tissue as a function of position, tiie fluorescence characteristic corresponding to at least one of 
fluorescence lifetime, fluorescence quantum efficiency, and fluorescence absorption, said 
processor being configured to generate an image signal as a fimction of said values; and 

(d) an output device responsive to said image signal to provide an image corresponding 
to the heterogeneous composition of the tissue. 

22. The system of claim 2 1 , fiirther comprising a number of modulated light sources. 

23 . The system of claim 2 1 wherein said fluorescence characteristic corresponds to 
fluorescent yield. 
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24. The system of claim 2 1 , wherein said sensor is configured to detect said emission at a 
number of locations along a surface of the tissue. 

25. The system of claim 21, wherein said processor is configured to determine said values 
fi-om a comparison of a calculated emission to an observed emission derived from said 
detected light signal, said calculated emission being determined as a function of an estimated 
spatial variation of said fluorescence characteristic, said estimated variation being updated 
and said comparison being repeated until a difference between said calculated emission and 
said observed emission reaches a desired minimum. 

26. The system of claim 21, wherein said source mcludes a laser diode and said sensor 
includes a CCD camera. 
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